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(57) Abstract 

A receiving apparatus Doppler spread estimator includes a receiver receiving signals over a transmission channel that varies over time 
due to fading and converting the signals to digital samples for processing. A stationarizer operatively associated with the receiver stationarizes 
the digital samples. An estimator is operatively associated with the stationarizer for estimating spectral density of the stationarized digital 
samples. The correlator is operatively coupled to the estimator for correlating the spectral density estimates with functions corresponding 
to hypotheses of plural Doppler spread values to select one of the Doppler spread values as an estimate of actual Doppler spread. 
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DOPPLER SPREAD ESTIMATION SYSTEM 

FIELD OF THE INVENTION 
This invention relates to wireless communication systems and, more 

particularly, to a method and apparatus for Doppler spread estimation in a wireless 

communication system. 

BACKGROUND OF THE INVENTION 

In wireless mobile communication systems the transmission channel can vary 
with time. This is due to the mobility of the transmitter and/or the receiver. Variations in the 
transmission channel are characterized by the Doppler spectrum. The Doppler spectrum is the 
Fourier transform of the autocorrelation function of the channel's time impulse response. The 
width of the Doppler spectrum is referred to as the Doppler spread in "Microwave Mobile 
Communications," W. C. Jakes, Jr., ed., New York: Wiley, 1974, IEEE Press, 1994, and 
"Wireless Communication," Theodore S. Rappaport, Prentice Hall PTR, 1996. The Doppler 
spread is commonly used as a measure of the channel variation rate or fading rate. 

It is often desirable to estimate the Doppler spread. For example, in order to 
adaptively optimize a coherent detector in a receiver, the parameters of the adaptation 
algorithm are made to be a function of the estimated Doppler spread. Also, in a cellular mobile 
telephone system, the handoff process can be enhanced if an estimate of the Doppler spread 
is available. This avoids handoff of fast moving users to micro cells. 

U.S. Patent Nos. 4,723,303 and 5,016,017 discuss the use of Doppler spread 
estimation. Both patents describe estimating the Doppler spread from the received signal itself 
However, these methods may not give a good estimate with a non-stationary received signal. 
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Lars Lindbom, "Adaptive Equalization for Fading Mobile Radio Channels," Techn.Lic.Thesis 
No. UPTEC 92124R, November 1992, Department of Technology, Uppsala University, 
Uppsala, Sweden, discloses a method of estimating the Doppler spread from a sequence of 
channel estimates. This method uses differentials of the channel estimate, which comprise 
5 difference of values between two points in time. However, the differential itself is generally 

very noisy and needs averaging. As a result, the average gives biased estimates of the Doppler 
spread. 

Still another proposal improves Doppler estimation by adaptive filters 
controlled by the final Doppler spread estimate. However, this is done with a significant 
10 increase of receiver complexity, while the estimate remains biased. 

The present invention is directed to overcoming one or more of the problems 
discussed above in a novel and simple manner. 



SUMMARY OF THE INVENTION 
In accordance with the invention there is provided a receiver Doppler spread 
1 5 estimator which utilizes stationarization of a non-stationary received signal. 

Broadly, in accordance with one aspect of the invention, there is disclosed a 
Doppler spread estimator including receive means for receiving signals over a transmission 
channel that varies over time due to fading, and converting the signals to digital samples for 
processing. First processing means are operatively associated with the receive means for 
20 stationarizing the digital samples. Second processing means are operatively associated with 

the first processing means for estimating spectral densities of the stationarized digital samples. 
Third processing means are operatively coupled to the second processing means for correlating 
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the spectral density estimates with functions corresponding to hypotheses of plural Doppler 
spread values to select one of the Doppler spread values as an estimate of actual Doppler 
spread. 

It is a feature of the invention that the receive means is adapted to receive code 
division multiple access (CDMA) signals. 

It is a further feature of the invention that channel coefficients are obtained 
directly from the CDMA signals. 

It is another feature of the invention that the receive means is adapted to receive 
signals including pilot symbols spread with a known CDMA spreading code. 

It is another feature of the invention that the received signal represents a 
sequence of transmitted symbols and the first processing means comprises a mixer for mixing 
the digital samples with a complex conjugate of the transmitted symbols. The complex 
conjugate of the transmitted symbols is derived from a channel estimator. 

It is another feature of the invention that the first processing means comprises 
a low pass filter. 

It is another feature of the invention that the second processing means develops 
vector data representing spectral components. 

It is still another feature of the invention that the third processing means 
comprises a multi-channel correlator. The correlator calculates correlations between estimated 
spectral densities and weighting functions corresponding to the hypotheses of the plural 
Doppler spread values. The third processing means includes a selector for selecting the 
hypothesis producing a highest correlation value. 
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In accordance with another aspect of the invention, there is disclosed a 
receiving apparatus Doppler spread estimator including a receiver receiving signals over a 
transmission channel that varies over time due to fading and converting the signals to digital 
samples for processing. A stationarizer operatively associated with the receiver stationarizes 
5 the digital samples. An estimator is operatively associated with the stationarizer for estimating 

spectral density of the stationarized digital samples. The correlator is operatively coupled to 
the estimator for correlating the spectral density estimates with functions corresponding to 
hypotheses of plural Doppler spread values to select one of the Doppler spread values as an 
estimate of actual Doppler spread. 

1 0 There is disclosed in accordance with a further aspect of the invention a method 

of estimating Doppler spread of a received signal comprising the steps of receiving the signals 
over a transmission channel that varies over time due to fading, and converting the signals to 
digital samples for processing; stationarizing the digital samples; estimating spectral density of 
the stationarized digital samples; and correlating the spectral density estimates with functions 

15 corresponding to hypotheses of plural Doppler spread values to select one of the Doppler 

spread values as an estimate of actual Doppler spread. 

More particularly, the present invention relates to a method and apparatus for 
Doppler spread estimation in a wireless communication system. The Doppler spread 
information can be used to optimize performance of a communication system. Applications 

20 include both base and mobile station receivers where a coherent detector is the preferred 

detector type. 

Further features and advantages of the invention will be readily apparent from 
the specification and from the drawings. 
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BRIEF DESCRIPTION OF THE DRAWINGS 
Figs. 1-3 are block diagrams of receiving apparatus including a Doppler spread 
estimator according to the invention; and 

Fig. 4 is a block diagram of the Doppler spread estimator according to the 

5 invention. 



DETAILED DESCRIPTION OF THE INVENTION 
Referring to Fig. 1, a receiving apparatus 10 according to the invention is 
illustrated. The receiving apparatus 10 is adapted for use with known pilot symbols. The 
receiving apparatus 10 includes an antenna 12 for receiving radio signals represented at 14. 

10 The antenna 12 is coupled to a radio receiver and converter 16. The radio receiver and 

converter 16 filters amplified signals and converts them to a suitable form for processing, such 
as complex numerical sample values. The signals are in the nature of digital samples for 
processing. The radio receiver and converter 16 is coupled to a channel estimator 18. The 
channel estimator 1 8 is connected to a Doppler spread estimator 20 and receives known 

15 symbols from a known symbol block 24. The known symbols comprise pilot or reference 

symbols also included in the received signal used for determining channel estimates. The 
known symbols are stored in a memory or generated locally in a code generator. The channel 
estimator 1 8 correlates the received digital samples with the known symbols to provide channel 
estimate to the Doppler spread estimator 20. The Doppler spread estimator 20 estimates the 

20 Doppler spread which is sent to a signal processing block 22. The signal processing block 22 

processes the sampled signal to extract information. The signal processing block 22 is also 
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connected to the channel estimator 1 8 so that after Doppler estimation the channel estimation 
can be improved. 

Referring to Fig. 2, a block diagram of a receiving apparatus 10' is illustrated. 
The receiving apparatus 10' is generally similar to the receiving apparatus 10 of Fig. 1. For 
5 simplicity, like elements are referenced with like numerals, while modified elements are 

indicated with primed numerals. Particularly, the receiving apparatus 10' differs in replacing 
the known symbols at block 24 with a symbol estimator 24'. This receiving apparatus is used 
in an application where symbols are unknown. The channel estimator 1 8 attempts to estimate 
the symbols, which are then stored in the symbol estimator 24'. The estimated symbols are 
10 subsequently used by the channel estimator 18 for determining channel estimates, as above. 

In the United States code-division multiple access (CDMA) cellular system 
known as IS-95, a transmitter transmits a stream of known symbols known as the pilot code. 
The pilot code is transmitted on the same channel at the same time as other, information 
bearing, symbols using different spreading codes. Fig. 3 illustrates a receiving apparatus 10" 
1 5 used in such a CDMA system. With CDMA systems the channel can be estimated directly to 

be used for Doppler spread. In this case, the known symbol block 24 of Fig. 1 is omitted. The 
channel estimator 18 correlates the received signal, comprising the pilot code and other codes 
in additive superimposition, with the known pilot code, and low pass filters the resultant 
complex correlation to obtain channel estimates. The received signal is also correlated with 
20 other codes carrying information it is desired to decode. The results of correlating with 

information carrying codes are multiplied by the conjugate of pilot code correlations for the 
same delay, and the results added to coherently combine the multi-path signals. 
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In wide band CDMA (WBCDMA) systems, modulation symbol intervals are 
much shorter. This allows multiple propagation paths to be resolved with much finer time 
resolution. 

In mobile communication systems, the received signal in a flat Rayleigh channel 
5 can be described by the following baseband model: 

r(t) = s(t) ■ c(t) + n(t), 

where 

s(t) = Ys&pft- n ^)> 

n 

d n is an n-th transmitted symbol, p(l) is a transmitted pulse, c(t) is a complex Gaussian process 
10 with mean zero and envelope correlation function 

ATr(r)=cr c 2 Jo(^r), 

&1 is a variance of the process c(l) % J^x) is a Bessel function of the first kind and of order 
zero, 

<o d = cOqv/Cq, co 0 is a carrier frequency, v is a mobile's velocity, c 0 is a light speed, and n(t) is 
1 5 white noise with the spectral density g. 

There are two conventional approaches for the basic Doppler spread estimation 
solution. These approaches are spectral analysis of the received signal or correlation analysis 
of the received signal. With spectral analysis of the received signal: 



1 Q ' x 



X r/exp(-y2^— ) 

i=qN iy 
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where 

n = dt- g + rii y 

r i9 c h and #i, are the i-th samples of the processes r(t) y c(t) y and n(t) y respectively. For this 
algorithm an estimate of the Doppler spread can be expressed in terms of FFT coefficients: 

r 

5 f*= -j^.arg[max5(*)] 

where/^ is a sampling frequency. 

With correlation analysis of the received signal: 

An estimate of the Doppler spread can be expressed in terms of the correlation function 
10 K(m) , calculated for some value m 0 : 

fo = k; x [k (mo)] 

where K^'Cc) is a function which is the inverse of the function K c (t). 

However, a vector r = (r„ r N ) is a time-varying, non-stationary random 
discrete process. Traditional spectral analysis does not work for non-stationary processes. To 
1 5 solve this problem, the present invention takes into consideration that d t = exp(J<f>^). Then, the 

algorithm 

fi - (d)* ■ n, 
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where d t is an estimate of the symbol d, , transforms the non-stationary vector r to the 

stationary vector f = (ri,...,w/) . 

Indeed, for non-dispersive channels dt » di . Besides, statistical characteristics 



of the noises n f and fit = di • n* are the same, because di 



2 



= 1 . Therefore, 



It is possible to improve the algorithm (9) by using the moving average: 



I iM+JVV-l 



n=iM 



In this case the errors Adn -d - dn are reduced and Equation (10) becomes more precise. 
For CDMA systems, channel coefficients (more exactly, channel coefficients estimates c i ) are 
10 directly obtained. In this case, channel coefficient estimates can be processed instead of /* . 

In accordance with the invention the non-stationary vector r is transformed to 
reduce the problem to analysis of the stationary vector f . Thereafter, the invention uses the 
optimal approach for estimating Doppler spread. This approach uses a criterion of the 
maximum of likelihood ratio. The optimal estimation algorithm is obtained by using the 
1 5 following equation: 
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^A(rl^) 1 



= 0, 



where A(r\ fd) is the likelihood ratio which is formed based on the vector f . 

Because the vector f has a Gaussian distribution, with the accuracy up to an 
inessential constant, the following equation can be written for the log-likelihood ratio: 



A(r 



k=0 



where G}f(k\fj) is a weight function which corresponds to some value of the Doppler 

spread f d . 

The method of estimating Doppler spread in accordance with the invention 
utilizes stationarization of the received signal, spectral analysis of the stationary vector, and 
correlation of the spectral estimate. Particularly, for stationarization of the received signal: 



Spectral analysis of the stationary vector f = (ri,...,rv> using a fast Fourier 
transform (FFT): 
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1 G-l 



(<7+l)"-l 

X * exp(-y2^-— )| 2 



Correlation is done in the frequency domain of the spectral estimate S(k) with 
L weight functions GH(k\fd) which correspond to L values of the Doppler spread f d \ 



A(r 



N/2-1 



Jt=0 



5 where f d = {fji\\ • is a vector of Doppler spread hypotheses. One of the Doppler 

spread hypothesis fj[m) is selected if: 

A(r|/*(m))> A(r\Mp)) 

Fig. 4 illustrates a block diagram of the Doppler spread estimator 20 of Fig. 1 . 
Stationarization of the received signal is performed using a multiplier 26 and low pass filter 28. 
10 The multiplier 26 receives the sampled received signals r n which are multiplied with the 

complex conjugate of the transmitted symbols dn . The radio receiver and converter block 

16, see Fig. 1, is the source of the received signal r„. The source of the complex conjugate of 
the transmitted symbols can be the channel estimator 1 8 or the known symbols 24. The output 
of the multiplier 26 is passed through a low pass filter 28 which reduces statistical errors 

1 5 Adn = dn — dn of the transmitted symbol estimation and out-of-band noise. 
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The output of the filter 28 is coupled to an input of a second processing block 
30. The block 30 estimates the spectral density of the incoming signal rt using the algorithm 
of Equation (15). The output, represented in bold lines, indicates that the output data is in the 
form of vector data which are the spectral components S = / 2)] . 

5 The output of the block 30 is fed to an input of a third processing function 32 

in the form of a multi-channel correlator which produces the likelihood ratio metrics using 
Equation (16). Each channel of the correlator calculates correlations in frequency domain 

between the spectral estimate S(k) and the waiting function Gx(k\fd). The waiting function 

Gh (k\fj(my) for the m-th channel corresponds to some selective value f/m) of the Doppler 

10 spread. Particularly, these comprise hypotheses of plural Doppler spread values. The 

likelihood ratio metrics are compared with each other using Equation (16). The outputs are 
applied to a max function block 34. On the basis of comparisons between the likelihood ratio 
metrics of L Doppler spread hypotheses, the multi-channel correlator 32 generates the decision 
that the Doppler spread is f/rn) if the m-th channel has the maximum output value. 

1 5 The illustrated method for estimating Doppler spread can be used in WCDMA 

or IS-95 systems, since in these systems known symbols or pilot symbols are periodically 
inserted into the data stream at every slot. 

As will be appreciated by one of ordinary skill in the art, the present invention 
may be embodied as methods or devices. Accordingly, the present invention may take the 

20 form of an entirely hardware embodiment, an entirely software embodiment, or an embodiment 

combining hardware and software aspects. The present invention has been described in part 
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with respect to the block diagram illustrations of Figs. 1-4. It will be understood that each 

block of the illustrations, and combinations of blocks, can be implemented by computer 

o 

program instructions. These program instructions, which represent steps, may be provided to 
a processor to produce a machine. 

5 Accordingly, blocks of the block diagrams support combinations of means for 

performing the specified functions in combinations of steps for performing the specified 
functions. It will be understood that each block of the illustrations, and combinations of 
blocks, can be implemented by special purpose hardware-based systems which perform the 
specified functions or steps, or combinations of special purpose hardware and computer 

10 instructions. 

Thus, in accordance with the invention, a Doppler spread estimator utilizes 
stationarization of non-stationary received signals and calculates correlation functions which 
are compared to hypotheses to determine an accurate and reliable estimate of Doppler spread. 
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WE CLAIM: 

1 . A Doppler spread estimator comprising: 

receive means for receiving signals over a transmission channel that varies over 
time due to fading and converting the signals to digital samples for processing; 
5 first processing means operatively associated with the receive means for 

estimating spectral densities of the digital samples; and 

second processing means operatively coupled to the first processing means for 
correlating the spectral density estimates with functions corresponding to hypotheses of plural 
Doppler spread values to select one of said Doppler spread values as an estimate of actual 
10 Doppler spread. 

2. The Doppler spread estimator of claim 1 wherein said receive means is 
adapted to receive Code-Division Multiple Access (CDMA) signals. 



3. The Doppler spread estimator of claim 2 wherein channel coefficients 
are directly obtained from the CDMA signals. 

1 5 4. The Doppler spread estimator of claim 1 wherein said receive means is 

adapted to receive signals including pilot symbols spread with a known CDMA spreading 
code. 
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5. The Doppler spread estimator of claim 1 further comprising third 
processing means operatively associated with the receive means for stationarizing the digital 
samples prior to estimating spectral densities of the digital samples. 

6. The Doppler spread estimator of claim 5 wherein said received signal 
5 represents a sequence of transmitted symbols and said third processing means comprises a 

multiplier for mixing the digital samples with a complex conjugate of the transmitted symbols. 

7. The Doppler spread estimator of claim 6 wherein said complex 
conjugate of the transmitted symbols is derived from a symbol estimator. 

8. The Doppler spread estimator of claim 5 wherein said third processing 
1 0 means further comprises a low pass filter. 

9. The Doppler spread estimator of claim 1 wherein said first processing 
means develops vector data representing spectral components. 

10. The Doppler spread estimator of claim 1 wherein said second processing 
means comprises a multichannel correlator. 
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11. The Doppler spread estimator of claim 10 wherein said correlator 
calculates correlations between estimated spectral densities and weighting functions 
corresponding to the hypotheses of the plural Doppler spread values. 

12. The Doppler spread estimator of claim 11 wherein the second 
5 processing means includes a selector for selecting the hypothesis producing a highest 

correlation value. 
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13. A receiving apparatus Doppler spread estimator comprising: 

a receiver receiving signals over a transmission channel that varies over time 
due to fading and converting the signals to digital samples for processing; 

an estimator operatively associated with the stationarizer for estimating spectral 
densities of the stationarized digital samples; and 

a correlator operatively coupled to the estimator for correlating the spectral 
density estimates with functions corresponding to hypotheses of plural Doppler spread values 
to select one of said Doppler spread values as an estimate of actual Doppler spread. 



14. The receiving apparatus Doppler spread estimator of claim 13 wherein 
said receiver is adapted to receive Code-Division Multiple Access (CDMA) signals. 



1 5 . The Doppler spread estimator of claim 1 4 wherein channel coefficients 
are directly obtained from the CDMA signals. 



1 6. The receiving apparatus Doppler spread estimator of claim 1 3 wherein 
said receiver is adapted to receive signals including pilot symbols spread with a known CDMA 
spreading code. 

17. The Doppler spread estimator of claim 13 further comprising a 
stationarizer operatively associated with the receiver for stationarizing the digital samples prior 
to estimating spectral densities of the digital samples. 
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1 8. The receiving apparatus Doppler spread estimator of claim 1 7 wherein 
said received signal represents a sequence of transmitted symbols and said stationarizer 
comprises a multiplier for mixing the digital samples with a complex conjugate of the 
transmitted symbols. 

19. The receiving apparatus Doppler spread estimator of claim 1 8 wherein 
said complex conjugate of the transmitted symbols is derived from a symbol estimator. 

20. The receiving apparatus Doppler spread estimator of claim 1 8 wherein 
said stationarizer further comprises a low pass filter. 

2 1 . The receiving apparatus Doppler spread estimator of claim 1 3 wherein 
said estimator develops vector data representing spectral components. 

22. The receiving apparatus Doppler spread estimator of claim 1 3 wherein 
said correlator comprises a multichannel correlator. 

23 . The receiving apparatus Doppler spread estimator of claim 22 wherein 
said correlator calculates correlations between estimated spectral densities and weighting 
functions corresponding to the hypotheses of the plural Doppler spread values. 
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24. The receiving apparatus Doppler spread estimator of claim 23 wherein 
the correlator includes a selector for selecting the hypothesis producing a highest correlation 
value. 
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25. The method of estimating Doppler spread of a received signal, 
comprising the steps of: 

receiving the signals over a transmission channel that varies over time due to 
fading and converting the signals to digital samples for processing; 
5 estimating spectral densities of the stationarized digital samples; and 

correlating the spectral density estimates with functions corresponding to 
hypotheses of plural Doppler spread values to select one of said Doppler spread values as an 
estimate of actual Doppler spread. 

26. The method of claim 25 wherein said receiving step receives Code 
1 0 Division Multiple Access signals. 



27. The Doppler spread estimator of claim 26 wherein channel coefficients 
are directly obtained from the CDMA signals. 

28. The method of claim 25 wherein said receiving step receives signals 
including pilot symbols spread with a known CDMA spreading code. 

15 29. The Doppler spread estimator of claim 25 further comprising third 

processing means operatively associated with the receive means for stationarizing the digital 
samples prior to estimating spectral densities of the digital samples. 
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30. The method of claim 25 wherein said received signal represents a 
sequence of transmitted symbols and said stationarizing step comprises multiplying the digital 
samples with a complex conjugate of the transmitted symbols. 

31. The method of claim 30 wherein said complex conjugate of the 
transmitted symbols is derived from a symbol estimator. 

32. The method of claim 30 wherein said stationarizing step comprises low 
pass filtering the mixed samples. 

33. The method of claim 25 wherein said estimating step develops vector 
data representing spectral components. 

34. The method of claim 25 wherein said correlating step comprises using 
a multichannel correlator. 

3 5 . The method of claim 34 wherein said correlator calculates correlations 
between estimated spectral densities and weighting functions corresponding to the hypotheses 
of the plural Doppler spread values. 

36. The method of claim 3 5 wherein the correlating step includes selecting 
the hypothesis producing a highest correlation value. 
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